Supersolid phase with cold polar molecules on a triangular lattice 
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We study a system of heteronuclear molecules on a triangular lattice and analyze the potential of 
this system for the experimental realization of a supersolid phase. The ground state phase diagram 
contains superfluid, solid and supersolid phases. At finite temperatures and strong interactions 
there is an additional emulsion region, in contrast to similar models with short-range interactions. 
We derive the maximal critical temperature Tc and the corresponding entropy S/N — 0.04(1) for 
supersolidity and find feasible experimental conditions for its realization. 

PACS numbers: 03.75.Hh, 67.85.-d, 64.70.Tg, 05.30.Jp 



Long-range interactions are a key ingredient in many 
models of strongly correlated electronic systems and frus- 
trated quantum magnets [l|. Magnetic dipolar interac- 
tions often occur in materials science and compete with 
short-range ferromagnetic interactions, which leads to 
spatially modulated phases The influence of long- 
range interactions is currently also attracting a lot of in- 
terest in cold atomic and molecular gases: the first sig- 
natures of long-range interactions have been observed for 
magnetic interactions in ^^Cr [3, 4], while electric dipole 
and van der Waals interactions between Rydberg states 
give rise to intriguing collective phenomena [5]. In addi- 
tion, there are big experimental efforts towards the real- 
ization of quantum degenerate polar molecules [sl, [t], [8| , 
where the permanent dipole moment of the molecules 
gives rise to a strong and hi ghly tunable electric dipole- 
dipole interaction [gj, [TO, [HI, 111 [11. 

In this Letter, we will concentrate on one intriguing as- 
pect of dipolar systems, namely the possibility to observe 
a supersolid phase in a single component system. While 
the interpretation of super-phenomena observed with tor- 
sional oscillators in solid ^He remains a puzzle [ll|, super- 
solids might be much easier to realize in lattices. While 
there exist a number of lattice models of hard-core bosons 
0, E 0, B 0,0,0, 0, soft-core bosons @ and 
quantum spins 0, [25[ most of these models are hard 
if not impossible to implement in a material. Systems 
of cold hetero-nuclear molecules are described by similar 
Hamiltonians but with longer range dipolar interactions. 
Here, we demonstrate that they show supersolidity under 
feasible experimental conditions. 

We consider bosonic polar molecules in a strong elec- 
tric field along the z-direction, which induces the dipole 
moment dz < d; here d denotes the permanent dipole 
moment of the hetronuclear molecule. The dominant in- 
teraction between the polar molecules is then given by 



the dipole-dipole interaction V{R) = 



the strength of the dipole interaction can be continuously 
tuned by the strength of the electric field. In addition, 
the polar molecules are confined into the xy plane by a 



where 
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FIG. 1: (Color online). Ground state phase diagram for the 
Hamiltionian Eq.([T]) around n = 1/3. The phases are a su- 
perfluid 'SF', supersolid 'SS' and a commensurate solid at 
density n = 1/3. With the double line we indicate a tran- 
sition region of the Spivak-Kivelson bubble type (emulsions) 
gradually going over to a region of incommensurate, floating 
solids with increasing interaction strength. For large interac- 
tion strength, and starting around half filling, the supersolid 
phase is suppressed by emerging solid ordering (stripes at half 
filling and incommensurate, floating solids at other fillings). 



strong transverse harmonic confinement as can be eas- 
ily achieved by a strong one-dimensional standing laser 
along the z direction. The combination of strong trans- 
verse tr app ing and dipole interaction creates a repulsive 
barrier ^| , which prevents the collapse naturally present 
in bosonic dipolar gases [26]. The effective long-range 
two-dimensional potential is then found by integrating 
over the z direction, and reduces to the effective 2D in- 
teraction V^^ D/r^. We refer to Refs. [TO, [H for a 
detailed discussion on how such a potential can be tai- 
lored. Subjecting the molecules to a triangular lattice 
created by three lasers in the xy plane, a standard one- 
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band tight binding analysis leads to an extended hard- 
core Bose-Hubbard model, 



-/i^n^. (1) 



Here, the first term describes the kinetic term with hop- 
ping amplitude /i is the chemical potential, ¥ = 0/0^ 
the 2D effective potential amplitude (a is the lattice spac- 
ing, and Ri denote the normalized lattice vectors). We 
will work with periodic boundary conditions, the hopping 
is set to unity, t = 1, the linear system size is L and the 
density denoted by n = N jl? . 

Our main results are the ground state phase diagram of 
Fig. [1] featuring a superfluid, commensurate solid and su- 
persolid phase and the finite temperature phase diagram 
for Vji < \h at constant density n = 0.4. The entropy 
SjN = 0.04(1) of the supersolid phase at the highest 
Tc is approximately one third of the entropy found in the 
system with short-range interactions at the same temper- 
ature and density, but is still comparable to the lowest 
entropies reached with bosonic ultracold alkali gases [271]. 
We find transitions that are continuous or belong to the 
Spivak-Kivelson bubble- type transition \2§\ . 

We study the Hamiltonian Eq. ([1]) by an unbiased 
and accurate quantum Monte Carlo simulation using the 
worm algorithm [29] in the implementation of Ref. [30]. 
To efficiently handle the slow (but integrable) decay of 
the potential in two dimensions, we replace the potential 
by a tabulated potential which is summed over all pe- 
riodic images (like in an Ewald summation) and which 
is only slightly different from In contrast to the 

short-range model, simulations with repulsive long-range 
interactions are notoriously more difficult and we suffer 
from the extra complication of a diverging number of low- 
energy metastable states [31]. Although this makes the 
identification of the phases and phase transitions at very 
large interactions V hard for current computational tech- 
niques, the supersolid phase can still be unambiguously 
identified in the experimentally most relevant regime. To 
identify the superfluid, solid and supersolid phases we 
measure the superfluid density ps and the density wave 

structure factor S^/L^ = {\Etli^ke'^^'^\^) / with 

Q = (4^/3,0). 

Ground state We start our analysis with the ground 
state phase diagram, shown in Fig. [1] for densities be- 
low half filling (the results above half filling are similar). 
This phase diagram has been obtained by extrapolating 
numerical data for different system sizes L = 12, 18, 24 
and sometimes L = 30 to infinity scaling the inverse tem- 
perature as /3t = L. The quality of the raw data from 
which the phase diagram was obtained can be assessed 
in Fig. O For V/t > 7.5(5) there is an insulating com- 
mensurate solid at density n = 1/3. 

For densities below 1/3 a superfluid phase is reached, 
similar to what is found for the short-range model [18, 19, 
[2Q| . [21I . I22I , but the tranisition here is different and of the 
bubble type introduced by Spivak and Kivelson [28] : over 




FIG. 2: (Color online). Density n, superfluid density ps and 
structure factor S^/L^ for F/t = 15 as a function of chemical 
potenial p. Statistical errors are smaller than symbol sizes if 
not shown. 
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FIG. 3: (Color online). Superfluid density ps and structure 
factor S^/L?' as a function of temperature T /tdX fixed density 
n = 0.4 and interaction strength V/t = 12. At low temper- 
atures we have a supersolid phase with both commensurate, 
solid order and a finite superfluid density. Although small 
system sizes show a continuous transition reminiscent of a 
two-dimensional 3-state Potts transition, a careful analysis 
reveals non-monotonicity and hysteresis in the structure fac- 
tor for larger system sizes, which is the first evidence for more 
complex physics in the emulsion phase. 
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a finite but narrow range of cilemical potentials, small 
crystallites form an emulsion of bubbles inside a liquid. 

For V = 30 (not shown) we find the first evidence 
for additional plateaus at various fillings below n = 1/3 
which are not present in the short-range model. With 
increasing system size the number of plateaus grows and 
they are separated by small superfluid regions. We ex- 
pect that an incommensurate, floating solid is formed in 
the thermodynamic limit for strong interactions by anal- 
ogy to the analysis of Ref. [32] . Note that in the classical 
limit of zero hopping the lon^-range model exhibits a 
devil's staircase (see Ref. jssi, [34| for Id) of various solid 
phases. 

Above the commensurate solid at n = 1/3 we find 
a continuous second-order phase transition belonging to 
the 3D XY model universality class to a supersolid phase, 
similar to what occurs in the short-range model. While 
near the tip the supersolid phase exists only over a narrow 
density range, it quickly extends {V/t = 15) all the way to 
half filling. For larger interactions {V/t > 20) and close 
to half filling, the structure factor and the superfluid 
density go down and supersolidity is lost for V/t = 30 at 
and near half filling. 

Finite temperature To study the the transitions at fi- 
nite temperature we will work in the canonical ensemble 
at a density n = 0.4. An interaction strength of at least 
V = 10.0(5) is needed in order to observe a supersolid 
phase (see Fig.|4]). For weak interactions we observe first 
a Kosterlitz-Thouless transition between a normal liq- 
uid and a superfluid phase, and a transition belonging 
to the two-dimensional 3-state Potts model then leads to 
the supersolid phase at lower temperature. When those 
two continuous transition lines cross for V/t = 11.8(5) at 
a temperature Tc = 0.53(8), the entropy per particle is 
S/N = 0.04(1) which is approximately one third of the 
entropy found in the short-range model with the same 
system parameters. For larger interaction strengths we 
observe in our simulations the emergence of an emulsion 
region with many metastable states between the normal 
liquid and the supersolid phase. More analysis for larger 
system sizes than what we can do in this study would 
be needed to accurately study the melting transition and 
the destruction of the superfluid order. 

Experimental proposal We now outline optimal pa- 
rameters for an experiment aiming at a homogeneous su- 
persolid phase in thermodynamic equilibrium. The ex- 
perimentally optimal density for the observation of the 
supersolid phase is given hy n ^ 0.4, as it exhibits the 
largest superfluid fraction and consequently the highest 
critical temperature Tc. The supersolid phase exists over 
a finite range of densities in the phase diagram, and this 
allows for flatter curvature in the trap center than with 
parabolic traps. 

From Fig. [H we see that the optimal interaction 
strength is V/t = 11.8(5) with Tc/t = 0.53(8). The re- 
pulsion should not be much larger than V/t = lb because 
of the risk of hitting a large emulsion region from which 
experiments will be unable to equilibrate due to the exis- 
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FIG. 4: (Color online). Finite temperature phase diagram 
at fixed density n = 0.4. The superfluid to normal liquid 
and the superfluid to supersolid transitions are found to be 
continuous (see text). After the Kosterlitz-Thouless (KT) 
and two-dimensional 3-state Potts transition lines cross, an 
intermediate phase of emulsions appears between the normal 
and the supersolid phases. The hatched region denotes where 
we found such emulsions for L = 24. 



tence of many metastable states. These parameters can 
be reached using LiCs with a dipole moment d = 6. 3D 
and an optical transition at A ^i:^ 940nm for the optical 
lattice [35]. Then, a polarization of the molecules with 
dz ^ O.ld gives rise to the interaction energy V ^ 0.23Er 
(here, Er ~ 1.6kHz denotes the recoil energy), and a rel- 
atively weak optical lattice with lattice ^ 8^r is sufficient 
to drive the system into the supersolid phase; fine tuning 
of the parameters can be achieved by controlling the in- 
teraction strength via the static electric field. Note, that 
only a weak polarization of the molecules is required, and 
consequently, the supersolid phase can also be reached 
for polar molecules with weaker dipole moments such as 
RbCs and LiRb. 

Finally, we note that standard time-of-flight images are 
an easy and direct tool to identify the different phases. 
In the superfluid phase the algebraic decay of the Greens 
function yields a strong signal at /c = and at all recip- 
rocal lattice vectors in the interference pattern. In the 
supersolid phase three times more peaks will show up 
since the unit cell contains three lattice sites. 

In conclusion, we have shown that fully polarized 
molecules loaded into a triangular optical lattice ex- 
hibit a supersolid phase. For weak interactions we have 
found in the ground state a bubble transition [28] be- 
tween a superfluid (n < 1/3) and a commensurate solid 
(n = 1/3) while a 3D XY transition leads to supersolid 
phase at higher densities when increasing the chemical 
potential. At finite temperature and fixed generic den- 
sity (say n = 0.4), there is a Kosterlitz-Thouless transi- 
tion from a normal liquid to a superfluid, and a 3-state 
Potts model transition to a supersolid. For larger inter- 
action strengths {V/t > 12.0(5)), there is a wide range 
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of temperatures in which an intermediate phase of emul- 
sions is found. The required temperature for reaching 
supersoUdity is feasible and the wide range in density 
allows for flexibility. We suggest an optimal value of 
V/t ^ 12 and a density of n = 0.4 when Tc ^ 0.53(8) or 
S/N = 0.04(1). Experimentally, the difference between 
the normal liquid, commensurate solid and supersolid 
phases can be detected by time-of-flight images. Our 
proposal is a very promising candidate for observing a 
clean supersolid phase in optical lattices. 
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